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Abstract 

We investigate Hamiltonian formulation of quantum scalar fields in a 
static quantum metric. We derive a functional integral formula for the 
propagator. We show that the quantum metric substantially changes the 
behaviour of the scalar propagator and the effective Yukawa potential. 

1 Introduction 

There is an old idea coming from Landau, Pauli and others ( see the review 
of the problem in 1957 by Deser [1]) that quantization of gravity could cure 
the ultraviolet divergencies in quantum field theory of matter fields interacting 
with gravity. Since then this inspiring idea has been followed by other authors 
(let us mention [2]- [3]). However, in general we find their approximate methods 
unreliable. In our earlier papers [4]- [6] we have investigated an interaction of a 
scalar field with a quantum or random metric in a functional integral formula- 
tion using the proper time representation for the propagator. The interaction of 
gravity with the scalar field has been treated non-perturbatively. We explicitly 
calculated an average over the metric field. For the reason of a mathematical 
consistency the functional measure of the metric field did not satisfy conditions 
required for the quantum field. In such a case it remained unclear whether the 
results could have been derived in a unitary framework of Hamiltonian dynam- 
ics. In this paper we study a similar model in the Hamiltonian formulation 
and relate the results to the ones of the functional integral form. We make an 
assumption that the metric does not depend on some coordinates in order to 
be able to perform explicit calculations. In particular, this assumption can be 
realized by a choice of a static metric. Then, we discuss in detail the model of 
a scalar field interacting with a static quantized metric. We are interested in 
a model where the quantum metric comes from a time-zero quantum field. In 
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general, such a field evolves with time. In order to simplify the model we ne- 
glect this time evolution. A perturbation of the free Hamiltonian by a quantized 
metric is a singular perturbation problem. So, wc first start with a regularized 
metric which is an analytic function. We discuss the operator Hamiltonian for- 
malism in the Fock space. We express the expectation values of the time-ordered 
products of the scalar field in the Fock space at some special complex points 
by means of a rigorous version of the Feynman integral. Then, the average 
over the metric field is calculated. Finally, the regularization of the metric field 
can be removed. We discuss in some detail an example of the metric being a 
threedimensional free quantum field. We show that an interaction with such a 
quantum metric substantially changes the short distance behaviour of the scalar 
propagator. The Hamiltonian approach allows to calculate the S-matrix and to 
calculate the effective two-body potential. We discuss these problems briefly. 



2 The Hamiltonian approach 

We consider the usual Lagrangian for the scalar field 

L = g^''dAct>dBcl> (1) 

Wc assume that the metric g in D dimensions (with the decomposition D = 
d+{D — d) of the coordinates x = (x, X))dcpcnds only on d spatial coordinates 
X and is of the block diagonal form (g^^) — ((/^'^(x), (/^"(x)) where g"^^ = 
rj^^ + h^^ and ry^^ = — 1 with the time t — Xd (we shall concentrate here on 
the case D = A and d = 2>). Let fa exp{—iuiat) be the set of positive frequency 
solutions (it is clear what it means because the metric is time-independent) of 
the wave equation 

-2Acl> = dAig^"" V9dB)4> = (2) 

where g = - det{gAB)- 

We choose the normalization 



J <^a 



We expand the solution </> into a complete set of positive fa exp{—ioJat) and 
negative energy solutions exp{ioJat) 

(j){t, x) = ^ /a(x)ac« exp(-ia;at) + /a(x)a+ exp(ia;ai) (3) 
where {aa, o+) are the usual annihilation and creation operators in a Fock space 

^ sc 

[aa-,a^]= 5al3 (4) 
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The Hamiltonian describing the time evolution (3) reads 

a 

Let |0)s be the Fock vacuum for the scalar field. We can compute the time- 
ordered product 

(0|T(<^(a:)<^(a:'))|0). = ^(/„(x)7„(x')+/c«(x')7jx)) exp(-za;„|t-i'|) ^ G^(x,a;') 

a 

(5) 

GiT is the solution of the equation 

dA{9''''V9dB)GF = -iS (6) 

It follows that {0\T{(f){x)(f){x'))\0)s = —iO~^(x,x'), where Dg denotes the wave 
operator on a manifold. As usual for the fields linear in creation and annihilation 
operators the higher order correlation functions are expressed by the two-point 
function. Then, the subsequent average over the metric concerns an average of 
the propagators Gp- 

We can also derive the time evolution (3) from the canonical formalism 

[</.(x),n(y)] =z5(x-y) 

and 

H = J dxVff(^(x)n(x)n(x)+.gJ"'^(x)a,(/)(x)afe(/)(x)) 

It is understood that there is a normal ordering of the bilinear forms in 11 and 
in H. The normal ordering requires metric-dependent counterterms. Such terms 
would change the gravitational dynamics but in our approximation wc have 
to neglect altogether the gravitational dynamics later on in order to perform 
explicit calculations. The non-linear local function of the quantum metric field 
in H itself needs a proper definition. The scalar field Hamiltonian leads to the 
wave equation (2) no matter what this definition is ( after the normal ordering 
of the scalar fields ). 

If g^^ = 1]^^ + h^^ then H = Hf) + Hj where Hj is the interaction Hamil- 
tonian for the computation of the .S'-matrix (see, eq.(8) below). 

The metric h^" may be operator-valued but we assume that the operators 
/i^^ commute (h^'^ are treated as a fixed classical background). In particular, 
we consider a model of relativistic fields h'^''{x,t) such that h^^iyi) are time- 
zero fields which live in the Fock space Tgr and are expanded in creation and 
annihilation operators ( we have in mind a model of a canonical free field) 

h'^^t, x) = ^ (/^"^(x) eM-ivgt)cg + exp(ii.„t)7r(x)4) (7) 
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where fa^i^) exp{—ii'at) are solutions of the wave equation on a certain classical 
background . The Hamiltonian for the metric field is 

a 

If the metric /; is defined on ii*^ (a flat background metric) then we have Ua = c|k| 
and a = {a, k) where cr is a polarization index, k is the momentum and 

/r = Pr (k)(27r)-^(c|k|)-i cos(kx) 

where P depends on the choice of coordinates (the gauge) [7]. We consider 

quantum gravitons at the temperature T described by the Gibbs density matrix 
exp{—(3Hg) where ^ = ksT and ks is the Boltzman constant. Then, the time 
ordered product of the metric fields reads 

Vi^^'^Pix- x') = (0|r(exp(-/3i?g)/i'"'(x)/iT''(x'))|0)g = 

{2Tr)-<^n J dk(c|k|)-ip^^;''^(k) cos(k(x-x'))exp(-ic|k||t-t'|)+ 

(27r)-''?i / (ik(c|k|)-ip'"';''''(k) cos(k(x - x'))(exp(c?i/3|k|) - ij-^ cos(c|k|(t - t')) 



where 



■ptJ-i^ilP — pMi' pjP 



depends on the gauge chosen for the metric field. In the limit h ^ (or 
T oo) we obtain the classical correlations resulting from the classical Gibbs 
distribution exp(— /3 / VhVh) of the initial values for the gravitational field. 
Hence, the covariance behaves as 



\x-x'\-c{t-t')\ 



-d+2 



With a general time- independent background metric we should replace c|k| by 
Vet and V cos(k(x — x')) by 

E(/r(x)7l''(x')+/j''(x')7r(x)) 

a 

in the formula for the temperature Green's functions . The short distance 
behaviour remains the same. 

Instead of the thermodynamic equilibrium described by the Gibbs distribu- 
tion we may consider a process of dissipation corresponding to an absorption of 
gravitational waves. Then, an absorption and emission of gravitational waves 
may lead to a time-independent (equilibrium) distribution of waves at large 
time, e.g., if the dissipation of gravitational energy is of the diffusive type, then 
for long times we could obtain 
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where the r-term describes a dissipation. If in the momentum (equiv. position) 
space the correlations of 6 have a 5 -type distribution then the invariant measure 
for the distribution of waves is the same as the one for the time-zero field at 
high temperature, i.e., A~^(x, y) ~ |x — y|~''+^. In any case in sec. 4 we assume 
the |x — yl""^"*"^ behaviour of the metric correlations no matter where it comes 
from. 

The S'-matrix can be calculated from its generating functional in the inter- 
action picture 

(0|T ( exp(-z J {Hi (s) + J dx J(s, x)</>(s, x))ds)) |0) 
= (0|det^-3exp(-iJGFJ)|0)g 

In this formula 11 and (p on the l.h.s. undergo the free Hamiltonian evolution 
but (as a simplification of our model) g does not change in time. In sec. 4 we 
calculate the expectation values of Gp under the assumption that g is Gaussian 
d = 3 and the propagator of g is |x — y|~^. From eq.(8) it can be seen that even 
if g is the free field then the expectation values over g of the Green's functions 
do not reduce to Gaussian integrals because of the g-dcpcndence of det A. In 
order to calculate an average over g we can apply the formula 

1 1 1 /" 

detA~i = exp(--Trln^) = exp(-Tr / dm^ {A + m^y^) 

subsequently expanding the exponential. Such an expansion corresponds to an 
expansion in the number of closed scalar loops. 



3 The Feynman integral for the propagator 

We apply a representation of the Feynman integral by means of stochastic pro- 
cesses [8] -[9]. For a general metric we would need a Brownian motion on a 
manifold for this purpose [10]. In order to simplify the argument we consider 
the wave operator on a manifold with a particular metric in D dimensions which 
in D = d + 1 = 4 can be related to the one called conformally static metric in 
refs.[ll]-[12]. We assume g^^^/g = 5^^ and goo = g" then 

-A = ^dAig^^'^dB) 

where g^^ = ^g^^ and h^^ = y/gh^^ , and we split the D-coordinates 
X = (x, X) into d coordinates x and D — d coordinates X. Eq.(2) in d = 3 can 
also be expressed in the form of a wave equation on a static thrccdimensional 
manifold (here A denotes the threedimensional Euclidean Laplacian) 

(^o' + ffooA).^ = (^o' + 5oo5fc(/'V590)'^ = 
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As an auxiliary tool for a calculation of the scalar propagator we consider 
the Schrodinger-type equation for the Hamiltonian A 

idrA (x, X) = Ai^r (x, X) (9) 

We can take the Fourier transform in X because g depends only on spatial 
coordinates x. Then, eq.(9) takes the form 

i^ri}{■^,P)=A^{■i^,P) (10) 

where 

A = -\/\ + \p,P.r''{^) (11) 

We consider the metrics h (as well as the initial states) which are analytic 
functions 

h{-x) = J dp/i(p) exp(ipx) (12) 

with the growth less than exp(e|zp) (with arbitrarily small e) for a complex z. 
Then, we can express the solution of eq.(9) by means of the Feynman integral 

^r{^,X)= j dPexp{iPX)E[Mri}{x + Xh{T),P)] (13) 

where 

X = Vi= ^(1 + i) 

and 

Mr = exp ( - '-P^P^ £ ff'"^(x + Xh{s))ds) (14) 
Eq.(14) is understood as a limit i? — > oo of a regularized expression with 

5fi(x + Ab(s)) = exp(-^)§(x + Ab(s)) (15) 

The proof of eq.(13) (for a Hamiltonian with a potential) was given in ref.[8]. 
The operator (11) coincides with the one for the Schrodinger equation with the 
potential \PgP. 

Eq.(13) can be considered as an analytic continuation of the imaginary time 
version (the diffusion equation [10]) 

-drM^.P)=Mr{^,P) (16) 

However, with = — 1 + wc would need to restrict ourselves to which 
are bounded from below if the diffusion equation (16) is to make sense. This is 
the basic reason for an analytic continuation from the Wiener to the Feynman 
integral in the final formulas in this paper. In contradistinction to the standard 
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methods we achieve the analytic continuation in time through an analytic con- 
tinuation in space. We would not need to require any analycity of the potential 
if we worked with the imaginary time. Wc would obtain the formulas (13) and 
(14) with A = 1. Then, an analytic continuation to complex A = "^(1 + *) is 
needed. Both the imaginary time and the real time formulas follow from the 
Trotter product formula . Denoting 

A = -\a + V (17) 

we have for the imaginary time 

(T r \ " 

exp(— A)exp(--y)j (18) 

The kernel K of the r.h.s. of the formula (18) can be expressed by the Brownian 
motion. In the real time case we express the kernel (for a finite n ) by means of 
the free propagators and rotate the integration axis from x to \x . Then, the 
limit n oo can again be expressed by the Brownian motion (as in eq.(13)). We 
can express the solution of eq.(13) by means of the Feynman kernel K defined 
by ^r{x) = J dyKr{x,y)tP{y). We have [8] [9] [13] 

Kr{x, y) = (27r)-^+'^(27riT)-i exp(^(y - x)^) / dPexp {iP {Y - X)) 
eM-iiP^)E[e^V {'Uo Pf^^''"' («' ^' y)) P-^s)] ^ ' 

where 

v(s,x,y) =x+-(y-x) + Ax/7a(-) (20) 
r T 

and a is the Brownian bridge starting from and ending at s = r in 0. The 
Brownian bridge is defined as the Gaussian process with mean equal to zero and 
the covariance 

E[aj{s)ak{s')] = 6jks{l-s') (21) 

for s < s'. The representation [13] a(s) = (1 — s)b(Y3j) (where b is the 
Brownian motion) is useful for computations. 



4 An average over the quantum metric 

We are interested in metrics which arc rclativistic quantum fields. In the real- 
time quantum field theory the time-ordered products have the Kallen-Lehmann 
representation [14] 

V^"'■'''f^{x,y) = {0\T{hf'^{x)h''l^{y))\0) = J dp(m^)Vp'''^ {x - y; m^) 
= Zime^o^"'''"^ / dp{m'^) idT{2mT)-i exp(-imV - er) exp(^(a; - yf) 

(22) 
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(23) 



Hence, 

= lirrie^o idTa{T) exp(-eT) exp(^(a; - yf) 
where 

cr(T) = J dp{m'^){2iTTT)~^ exp{-^m^T) (24) 
The Euchdean version of eq.(23) reads 

P(x - y) = XT drai-ir) exp(-^(x - y)^) (25) 
It can be seen that if 

x-y ^x-y+z 
then we can continue to complex z if 

K((x-y + z)2)>0 

We can formulate the analycity also in terms of the radial analycity meaning 
that X— y A(x— y) if 3(A-^(x— y)^) > 0. From the point of view of the proper- 
time representation in the imaginary time formulation it is useful to introduce 
the regularized imaginary time metric fields in the form 



poo ^ 

hs{x) = J rfrexp(--k^)x(k,r)-v/CTPir)exp(ikx) 



(26) 



where 



(F"(k, r)x''"(k', r')) = V^'''P''5{k + k')5(r - r') (27) 

Then, hs{x) is an analytic function of x and its correlation function is Vg where 

in cq.(25) the T-intcgral starts from t > S. For further computations wc assume 
that the metric is Gaussian. We have seen in eq.(8) that an average over the 
scalar field leads to a determinant of A which depends on the metric g. It 
is difficult to calculate non-perturbativcly a contribution of the non-Gaussian 
terms in the complete theory of the interaction of gravity with a scalar field. 
We are able to calculate the averages over the metric only through an expansion 
in the number of closed scalar loops as discussed at the end of sec. 2. In such an 
expansion the ultraviolet behaviour does not depend on the number of loops. 
Hence, we may restrict ourselves here to the zeroth order corresponding to 
det^ = 1. 

Using the Trotter product formula (18) (with V = \PgP and r —ir) we 
obtain (for a regular Gaussian metric (26), see eq.(8) for a calculation of the 
Gaussian average) 

(0| ( exp(i^ A) exp(-^y)) "|0),,(X, x; Y, y) 

= /dPexp(zP(X - Y)){2ml)-i exp(m^^^^) (28) 

exp - ^Y.PPVs{xj - Xj_i)PP^dxi....dx„ 
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In this equation we rotate the integration line x Ax. Then, the limit n oo 
can be expressed by the Brownian motion. As applied to propagators instead 

of the kernel A~^{x, y) consider ^"'^(Ax, X; Ay, Y). Using the Brownian bridge 
representation we obtain for an average (28) over the metric (22) 

(^-i(Ax, X: Ay, Y)) = J^"" dr cxp(zP(X - y))(2^T)-^ exp(-^(x - y)^) 
£;[exp ( - i /; /; dsds'PPVs{Xms) - Hs'))Pp)] 

(29) 

where 

s s 
v(s) = x + -(y - x) + v^a(-) 

T T 

We can calculate higher order correlations 

(>l-HAx, X; Ay, y)^-i(Ax', X'; Ay', Y')) = dn dT2dPdP' 
cxp(iP(X - Y))i2nn)-i exp(-2i^(x- y)2)exp(iP'(X' - Y')){2wT2)-i 
exp(-^(x'-y')2) 

£:[exp ( - i J^' J^' dsds'PPVs{X{^{s) - v{s'))PP 
-IJT J? dsds'P'P'Vs{X{^'{s) - v'{s'))P'P' 
+ i j;^ dsds'PPVs{X{ir'{s) - v{s'))P'P')] 
+{x x') 

(30) 

Here, {x x') means the same expression with x exchanged with x' , v' is an 
independent bridge with (x, y) — > (x',y'). In ref. [15] we have shown that the 
formulas (29)- (30) come from a resummation of the perturbation series in h of 
□ It follows from eq.(8) that this is also a resummation of the Dyson series 
in the interaction Hamiltonian. We could derive the short distance behaviour 
of the 4-point function from eq.(30). It can be shown [5] that its singularity is 
a product of the singularities of the two-point functions. 

Whether the correlations (29)- (30) are finite or not depends on whether the 
argument in the exponent is a well-defined function and whether the exponential 
factors are integrable. The functions in the exponent are analytic continuations 
from real arguments. It will be useful to use a concrete representation of these 
functions. So, for our special case of the |x— y|~^ metric correlations we consider 
the representation 

Then, we can perform the integral over k (with the regularization (26)) 

I?5(x - y) = /,°° dr / dkexp(-ir|k|2) exp(ik(x - y)) = 

/;° dr{2T,r)-i exp ( - \r-\^ - yf^ ^^'^> 
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From the representation (32) it can be seen that this is an analytic function of 
= (x — y)^ as long as 5fz^ > (the equality sign is allowed only if Sjz^ ^ 0). 

The argument of the exponential factors in cqs.(29)-(30) is a well-defined random 
variable because it is square intcgrablc. In fact, the expectation value 

E[\ J dsds'vs{\i ^s)-^s'm 

= fdsds'dtdt'E[Vs{X{^{s) - w{s'))Vs{X{w{t) - v(t'))] ^ ' 

can be explicitly calculated (using eq.(21)) and shown to be finite for (5 > as 
well as in the limit 5 — > 0. When S = the complex argument in the exponent 
in eqs.(29)-(30)) for the propagator (32) scales as 

-D(A(v(s) - v(/))) = -jV{v{s) - v(s')) (34) 

where T> is positive. Hence, the real part of (34) is negative supplying a damping 
factor for the integrals of the exponential factors (29)-(30). 

After the proof that the expressions (29)- (30) are finite we obtain just by 
scaling ( as in [5], although g^" = r]^" + hi^" is not scale invariant) 

(^-i(X,y)) ~ (35) 

at short distances (and in the approximation deiA = 1). 

For the behaviour in the x direction we put X = Y = Q then 

(^-i(x,y))=ii|x-y|-2+i (36) 

where is a constant. It is more regular than the two-point function for the 
four-dimensional free field (see ref.[16] for a more general argument). In order 
to treat the general case of 2n-point functions and det Iwe expand det A 
in eq.(8) in terms of the Green's functions as discussed at the end of sec. 2. 
Then, at each order we can calculate the expectation values of products of 
Green's functions (as in eq.(30), see also [5]). In each order of the expansion of 
det^l the results (35)-(36) can be generalized to arbitrary 2n-point correlation 
functions proving their anomalous short distance behaviour. 



5 Discussion 

We have discussed a Hamiltonian formulation of an interaction of a quantum 
metric with a quantum scalar field. We have compared the standard Dyson 
expansion for this model with the results coming from the functional integral. 
We have calculated scalar field correlation functions in the Gaussian approxi- 
mation for the metric field expanding the correlation functions in the number of 
closed scalar loops. The ultraviolet behaviour does not depend on the number 
of loops. We have shown that at each order of the expansion the short distance 
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behaviour is substantially modified as a result of the interaction with gravitons. 
In particular, the behaviour at equal times of the correlation functions is more 
regular than the one for the free scalar field. In sec. 4 we concentrated on the 
static metric corresponding toD = d+ l = 4of the framework of refs.[4]-[6] 
because it may be relevant to physical models of gravitons in equilibrium with 
matter. The Hamiltonian framework of this paper may be considered for any 
d < D. From the S-matrix formula we can obtain (inserting static sources in 
eq.(8); see also a more detailed discussion in [14] ) the expression for the Yukawa 
potential 

/oo 
dXD{A-\x,X,XD)) 
-oo 

here we denoted X = (X, Xn), where Xd is the time. If ©(x) ~ |x|~^''' at short 
distances then using just scaling properties of the formula (29) (as in refs.[4]-[6]) 
we can conclude that 

y(x,x = o)^ |x|2-<^-(i-^)(^-'^-i) 

Hence, in general the Yukawa potential is less singular in the x direction than the 
canonical one ( except the static case of sec. 4 corresponding to D = d+1 when 
the Yukawa potential does not change in the x direction at short distances). 
In the X direction {if D > d+1) 

y(x = o,x) = |x|-^+t^-('^+i)t^ 

The distance scale at which the anomalous behaviour would appear is deter- 
mined by the scaling behaviour of the metric field. In Einstein gravity there is 
a dimensional parameter of the Planck scale. In this case the eventual change 
of the short distance behaviour is expected below this length scale. 
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